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ABSTRACT 



Or 



Context. The determination of interstellar abundances is essential for a better understanding of stellar nucleosynthesis and the "chem- 
ical" evolution of the Galaxy. 

Aims. The aim is to determine l8 0/ 17 abundance ratios across the entire Galaxy. These provide a measure of the amount of enrich- 
ment by high-mass versus intermediate-mass stars. 

Methods. Such ratios, derived from the C 18 and C 17 7=1-0 lines alone, may be affected by systematic errors. Therefore, the C ls O 
and C 17 (1-0), (2-1), and (3-2), as well as the 13 CO(1-0) and (2-1) lines, were observed towards 18 prominent galactic targets (a 
total of 25 positions). The combined dataset was analysed with a large velocity gradient model, accounting for optical depth effects. 
Results. The data cover galactocentric radii between 0. 1 and 16.9 kpc (solar circle at 8.5 kpc). Near the centre of the Galaxy, 18 0/ 17 
= 2.88+0. 1 1 . For the galactic disc out to a galactocentric distance of ~ 10 kpc, 18 0/ I7 = 4.16+0.09. At ~ 16.5 kpc from the galactic 
centre, 18 0/ 17 = 5.03+0.46. Assuming that ls O is synthesised predominantly in high-mass stars (M>8M Q ), while C I7 is mainly a 
product of lower mass stars, the ratio from the inner Galaxy indicates a dominance of CNO-hydrogen burning products that is also 
' apparent in the carbon and nitrogen isotope ratios. The high 18 0/ 17 value of the solar system (5.5) relative to that of the ambient 

interstellar medium suggests contamination by nearby high-mass stars during its formation. The outer Galaxy poses a fundamental 
problem. High values in the metal-poor environment of the outer Galaxy are not matched by the low values observed towards the even 
more metal-poor Large Magellanic Cloud. Apparently, the outer Galaxy cannot be considered as an intermediate environment between 
the solar neighbourhood and the interstellar medium of small metal-poor galaxies. The apparent 18 0/ 17 gradient along the galactic 
disc and the discrepancy between outer disc and LMC isotope ratios may be explained by different ages of the respective stellar 
populations. More data from the central and far outer parts of the Galaxy are, however, needed to improve the statistical significance 
' of our results. 
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<^> 1. Introduction nuclear reaction rates, these reaction rates remain uncertain (see 

e.g. Stoesz & Herwig 2003). To summarise, measured 18 0/ 17 
Isotope abundance ratios play an important role in our under- ratios have the pote ntial to provide relevant constraints for mod- 
standing of stellar nucleosynthesis and the secular "chemical" els of stellar nuc i e osynthesis and "chemical" evolution, 
i enrichment of the interstellar medium by stellar ejecta. While 
Cd isotope ratios are not easily measured at optical wavelengths, ob- Th e 18 0/ 17 ratio is readily determined from the C 18 0/C 17 
servations of molecular clouds with their large number of molec- line intensity ratio. Both C ls O and C n O have similar chemical 
ular species allow us to distinguish between various isotopic and excitation properties. Rotational lines of both isotopologues 
species. In principle, accurate line intensity- and abundance- are traditionally assumed to be optically thin, although C 18 can 
ratios can then be determined. reach considerable optical depths in some regions (e.g., Bensch 
A particularly useful tracer of nuclear processing and metal et aL JHSJ- Oxygen isotopes are not fractionated (e.g., Langer 
enrichment is the I8 0/ 17 ratio because I8 must be released et aLH2S9, which is related to the high first ionization potential 
primarily from high-mass stars (M>8M ), while I7 may pre- of ox yg en ' 13 6 eV Thus ' ° + should have a low abundance in 
dominantly be ejected from stars with lower mass (Henkel & molecular clouds and charge exchange reactions like those for 
Mauersberger 1993; Henkel etal. 1994; Langer & Henkel 1995 ). carbon < see Ec l- 1 ln Watson et d-MIB should not significantly 
Heger & Langer (2000) concluded that at the surface of rotating affect atomlc and molecular oxygen abundance ratios. To sum- 
stars with masses of 8-25M Q , n O is enriched and 18 is depleted mans 7 e ' observed CO line intensity ratios are a good measure of 
during hydrostatic burning prior to the supernova event. Most ^/ O. 

mass is ejected later, during the supernova explosion. While Ratios of 18 0/ 17 have commonly been determined for 

Hoffman et al. (120011 ) conclude that production of both ls O and clouds in the solar neighbourh ood an d the inner G alaxy (see the 

17 in massive stars is reduced when taking into account new reviews of Wils on & Ma tteucc i [l992l H enkel e t al. [l994l Wilson 

& Rood [19941 Kahane [19931 Bieging [19971 ). For the galactic 

Send offprint requests to: J.G.A. Wouterloot e-mail: disc and centre region, Penzias dl981| ) reported average 18 0/ 17 

j .wouterloot@jach.hawaii.edu ratios of 3.65±0.15 and 3.5+0.2, respectively. He found no sig- 

* Figure fAjJ (the spectra) is available in electronic form nificant gradient with galactocentric distance out to kpc 
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(# =8.5 kpc). The 18 0/ I7 ratio from HCO + 7=1-0 in SgrB2 
(3.1+0.6; Guelin et al.[1982) is within the errors consistent with 
the CO results for the galactic centre region. The observed inter- 
stellar ratios as well as many ratios from envelopes of late-type 
stars of intermediate mass (some of them with 18 0/ 17 0<1; e.g. 
Kahane et al. 1992) differ significantly from the solar system 
18 0/ 17 ratio of 5.5, and Prantzos et al. ( 1 19961 ) even considered 
the possibility of measurement errors in the galactic data. 

The above-mentioned C ls O and C n O observations were 
made in the 7=1-0 transitions alone, which is not sufficient when 
trying to account for radiative transfer effects. There exist also 
some 7=2-1 data (e.g. Hofner et al. 120001 White & Sandell 
1 19958 ), but these have been used to analyse C 18 optical depths 
and column densities, rather than studying the isotope ratio it- 
self. Bensch et al. (2001) observed the definitely optically thin 
13 C ls O and 13 C n O 7=1-0 lines and derived 4.15+0.52 at one 
position in the nearby (c/~140pc) p Oph cloud. The same cloud 
was observed towards 21 positions in the 7=1-0, 2-1, and 3-2 
transitions of C I8 and C 17 by Wouterloot et al. (120051 1. Large 
Velocity Gradient (LVG) model calculations combining these 
data resulted in a similar but more accurate value, 4.11+0.14, 
for the I8 0/ I7 ratio. Very recently, Zhang et al. (2007) mea- 
sured the 7=1-0 lines towards a selected area of the southern 
star forming region NGC6334 (c/~1.7kpc). Under the (in this 
case realistic) assumption that all C I8 lines were optically thin, 
they obtained a ratio of 4.13+0.13. 

To date, little has been done to analyse molecular abun- 
dances in clouds at the edge of the Galaxy (for a pioneering 
study, see Wouterloot & Brand 1996). With respect to the so- 
lar neighbourhood, the far-outer Galaxy, at galactocentric radii 
of /?gc>16 kpc, is characterised by a low metallicity, a weak in- 
terstellar radiation field, a small cosmic ray flux, and a low large- 
scale average gas pressure (Brand & Wouterloot 1995 Ruffle et 
al. 120071 and references therein). 

Our initial measurements of two objects at large galactocen- 
tric ra dii, DDT94 Cloud 1 and DDT94 Cloud 2 (Digel et al. 
1199411, w ith distances of R G c of 22 and 17+2 kpc (Smartt et 
al. [1996 ) did not provide conclusive results. Therefore, to de- 
rive more reliable values, we have observed four objects from 
the catalogue of Wouterloot & Brand ( 1989) with stronger C I8 
lines (7^ ~ 1 K versus 0.2 K for the DDT94 clouds) to cover the 
range of Rgc beyond 15 kpc. In addition, we have reobserved 
objects at kpc ;£ Rgc ^10 kpc to compare these results with 
previous estimates based on the 7=1-0 line alone. Our sample 
consists of 25 positions in 18 sources that are observed in the 7 
= 1-0, 2-1, and 3-2 rotational lines of C ls O and C I7 0. I3 CO(l- 
0) and (2-1) data were also taken. 

2. Observations 

2.1. IRAM30-m 

On March 4, 1995, we used the IRAM 30-m telescopsfl to ob- 
serve C ls O and C I7 (1-0, 2-1) towards DDT94 Cloud 1 and 2. 
Between July 28 and August 2, 1995, we used the same tele- 
scope to observe C ls O (1-0, 2-1) and C 17 (1-0, 2-1) (see 
Table [TJ towards the other clouds of our sample. We used three 
SIS receivers simultaneously (first for C ls O 7=1-0, 2-1 and 
C n O 2-1, then for C n O 1-0, 2-1 and C ls O 2-1; all data were 
consistent) in combination with an autocorrelator split into three 
parts with equal velocity resolutions (in most cases 0. 10 km s _1 ). 



Table 1. Observed transitions. 
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12 C 17 O(l-0) 
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12 C ls O(3-2) 


329330.545 


14 
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12 C n O(3-2) 


337061.129 


14 
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" : Half Power Beam Width 

b : Beam (B cff ) and forward hemisphere (F c(f ) efficiency. To convert 

antenna temperatures to r mb main beam brightness temperatures, 
multiply the 7=1-0 and 2-1 values by F^/B^g (see, e.g., 
Rohlfs & Wilson 1996). For the 7=3-2 lines, multiply by l/5 eff . 



Towards most positions we initially used frequency-switching 
(by 7.9 and 15.8 MHz at 3 mm and 1.3 mm, respectively). 
Towards sources with broad lines, however, we observed in a 
position-switching (total power) mode. Offsets of the reference 
positions are given in Table [2] Velocity resolutions and rms 
values of the spectra are displayed in Table 3 (see Sect. 3). 
At positions with weak C I7 lines, frequency-switched spec- 
tra showed prohibitively large baseline ripples. Therefore these 
sources were also observed in total power mode. From most 
spectra obtained in total power mode we subtracted a low or- 
der baseline. Remaining baseline ripples due to standing waves 
were removed by subtracting a sinusoidal baseline. 

To investigate the importance of radiative transfer effects we 
observed 13 CO (1-0) and (2-1) simultaneously with C ls O (1-0) 
and (2-1) between May 10 and 15 and on September 7, 2000. 
All observations were made in a total power mode. For the four 
sources from the sample of WB89, the I3 CO and C I8 data were 
taken from Wouterloot & Brand ( 1996 ). 

On April 4-5, 2005, we used HERA to observe C I8 (2-1) 
towards all sources. HERA consists of two 9 pixel receivers on 
3x3 arrays (a 24" raster). The spectral resolution was 0.1 1 km/s 
and observations were made using position switching. These 
data are used to convolve the 7=2-1 and 7=3-2 spectra to the 
IRAM 7=1-0 angular resolution. 

An rms pointing accuracy of about 3" was determined from 
scans across nearby continuum sources. Planet measurements 
showed that the receiver alignment was within 1^1". The ob- 
served transitions, frequencies, beamsizes, beam and forward ef- 
ficiencies are given in Table Q] In the following, T* A intensities 
are given, unless specified otherwise. 

2.2. JCMT 15-m 



1 IRAM is supported by INSU/CNRS (France), the MPG (Germany), 
and the IGN (Spain). 



On July 14, 2001, we observed C n O(3-2) towards W49N 
and W51d with the 15-m James Clerk Maxwell Telescope 



Wouterloot et al.: Interstellar ls O/ 17 ratios 



(JCMTj^ using position switching. On various occasions be- 
tween November 7, 2004, and May 19, 2005, we observed 
C I7 0(3-2) towards NGC2024, L134N, and four offset positions 
near Ori-KL. C ls O(3-2) was observed towards the same sources 
and towards WB89-501. On February 10, 2006, we observed the 
remaining sources (except for the WB89 objects and G34.3+0.2, 
W30H, and NGC7538) in C ls O and C 17 0(3-2). 

The JCMT observations were made with the 2x800 chan- 
nel DAS autocorrelator and the dual polarization RxB3 SIS re- 
ceiver using position switching. The channel spacing was 156 - 
625 kHz and the bandwidth was 250 - 920 MHz, depending on 
the line width of the source. Offsets of the reference positions 
used are given in Table [2] Pointing was done on nearby contin- 
uum sources and was generally accurate to within 3". 

From the JCMT archives we obtained C 17 0(3-2) and 
C ls O(3-2) data towards G34.3+0.2, W30H, andNGC7538. The 
JCMT archival data are from standard sources which are fre- 
quently observed, and always after doing a pointing observa- 
tion on the source itself. The coordinates are within a few arcsec 
(much less than the beamsize and within the pointing accuracy) 
of the positions used for our IRAM 30-m observations. In some 
cases several spectra were available. Discrepant results and mea- 
surements with too low velocity resolution, affecting peak inten- 
sities, were not considered and the remaining spectra were aver- 
aged. 

2.3. Calibration uncertainties 

Calibration uncertainties in individual spectra may amount to 
+ 10% which corresponds to +14% in line intensity ratios. 
Changes in telescope efficiencies (see Table 1) may also con- 
tribute to the error budget. Like the uncertainties in the effi- 
ciencies, pointing errors may also affect those intensity ratios 
where the two lines were measured separately. This holds for all 
JCMT observations. As mentioned in Sect. 2.2, however, results 
from repeatedly observed sources turn out to be consistent. A 
direct way to assess the quality of the IRAM data is provided 
by Table 3. Here, line parameters for the C ls O 7=1-0 and 2-1 
lines are given twice, because the lines were once observed to- 
gether with C 17 and another time with 13 CO. The uncertainty 
in the IRAM line ratios will be less than those in the line inten- 
sities which are suggested by differences between both sets of 
C I8 measurements, because the ratios are commonly based on 
simultaneous observations. Nevertheless, the data displayed in 
Table 3 provide firm upper limits to the full error budget of the 
IRAM line intensities. 



3. Results 

Table [2] contains the sample of observed sources. SgrB2 is the 
most prominent star-forming region close to the galactic centre, 
L 134N is a dark cloud near the Sun, and WB89 380, WB89 391, 
WB89 437, and WB89 501 are sources located in the far outer 
Galaxy. The other targets are star forming regions located at an 
intermediate range of galactocentric radii, extending from the 
inner "molecular ring" out to the Perseus arm. 

2 The JCMT is operated by The Joint Astronomy Centre on behalf of 
the Science and Technology Facilities Council of the United Kingdom, 
the Netherlands Organisation for Scientific Research, and the National 
Research Council of Canada. Our programs: m04bdl0 and m05bil7. 
JCMT archive data were obtained from the Canadian Astronomy Data 
Centre, which is operated by the Herzberg Institute of Astrophysics, 
National Research Council of Canada. 
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The complete set of spectra is displayed in the appendix 
(Fig. I A. It that is available in electronic form. All spectra of an 
individual source are plotted on the same velocity scale. Line pa- 
rameters obtained from Gaussian fits to the 7=1-0, 2-1, and 3-2 
transitions are summarised in Table [3] To compare abundances 
and not line intensities, we corrected the C 18 0/C I7 line inten- 
sity ratios for the difference in frequency of both isotopologues, 
following Linke et al. ( 1 19771 1 and Penzias (1981). This implies 
an increase of the ratios by a factor 1.047. For W51d, all param- 
eters were calculated for two velocity intervals, accounting for 
the two emission lines at +50 and +60kms~ 1 detected towards 
this source. The new C ls O HERA line parameters towards the 
WB89 sources agree with those obtained earlier, also with the 
30-m telescope, by Wouterloot & Brand dl996l l. 

The resulting C I8 0/C 17 abundance ratios as a function of 
distance from the galactic centre, Rqc are shown in Fig. Q] 
Error bars accounting for the uncertainty in the line areas (ob- 
tained from baseline fits and given in Table [3), are not plotted to 
avoid confusion. For comparison we also show the results from 
Penzias (1981) for the 7=1-0 transition. For the 7=2-1 lines, 
we plot results from Hofner et al. (120001 ). who used the 30-m 
telescope to observe C ls O and C 17 towards Hn regions of the 
inner Galaxy as well as towards W3(OH). The good agreement 
between our results and those of Penzias (1981) and Hofner et 
al. J2000 ) indicates that the calibration of our data is correct. 

The distance to a given source can lead to a number of ob- 
servational biases that are difficult to quantify and that do not 
facilitate a comparison between different parts of the Galaxy. 
With increasing distance, the linear beam size increases, leading 
potentially to the inclusion of a significant amount of relatively 
diffuse optically thin gas. On the other hand, a tendency to in- 
clude brighter and more massive sources at larger distances may 
involve sources with systematically higher C I8 opacities that 
may lead to underestimates of the C 18 0/C I7 abundance ratios. 
To search for a distance-related bias, Fig. |2]therefore shows the 
C 18 0/C I7 ratios as a function of solar distance d. In Fig. [3] we 
show the same as in Fig.Q] but for the 13 CO to C I8 ratios ob- 
tained from the 7=1-0 and 2-1 lines. In addition, Figs. [3}; and 
d show correlations between the I3 CO/C 18 and C 18 0/C 17 ra- 
tios. 



4. Discussion 

4.1. Overall trends 

Not accounting for optical depth effects and omitting Sgr B2 
and the far outer Galaxy clouds, our average C 18 0/C 17 abun- 
dance ratios are 3.76±0.16 (7=1-0), 3.65+0.16 (7=2-1), and 
3.08+0.12 (7=3-2), which can be compared with the Penzias 
([198TT ) 7=1-0 galactic disc value of 3.65±0.15. The 7=1-0 
results agree within the errors. If we weight our data with 
their uncertainties, the average ratios become slightly smaller, 
i.e. 3.68+0.11, 3.38+0.13, and 2.74+0.10, respectively. We 
note, however, that smaller values tend to have smaller errors. 
Therefore, weighted averages may yield too small values and 
unweighted averages should be preferred. For the two SgrB2 
positions average ratios are 2.99+0.06 (1-0), 2.26+0.17 (2- 
1), and 2.47+0.38 (3-2) (weighted 2.99+0.05, 2.27+0.14, and 
2.32+0.30). The average ratios of the four far outer Galaxy 
clouds are 5.18+0.96 (1-0) and 5.02+0.22 (2-1) (weighted: 
4.79+0.77 and 4.96+0.23). 

Assuming that the C ls O and C n O rotational lines are opti- 
cally thin, the most important results that can be derived from 
Figs. Q] and [2] are: (1) Our C l& 0/C u O abundance ratios of the 
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Table 2. Observed sources". 
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(1) Brand & Wouterloot d 1994k peak position from Wouterloot & Brand d!996t . 

(2) Kinematic distance from the observed radial velocity and the rotation curve of Brand & Blitz ( 1993 1. For W33 and G34. 3+0.2, the near 
kinematic distance is given. 

(3) Weighted average of distances of Hachisuka et al. (2006), 2.04±0.07 kpc, Xu et al. (2006), 1.95+0.04 kpc, and Imai et al. (2000), 
1.83±0.14kpc. 

(4) Smartt et al. ( fl996l 

(5) Menten et al. d2007b. 

(6) Anthon y-Twar og (1982}. 

(7) Franco (1989). 

(8) Kerr & Lynden-Bell (1986 1; we note that Eisenhauer et al. (2005 1 derived 7.62+0.32 kpc for the distance to the Galactic Centre which may be 
compared with the IAU value of 8.5 kpc assumed in this work. 

(9) Gwinn et al. (l99"2l 

(10) Lacy et al. (2007) 

(11) Dickel et al. (1978}. DR 21cal position from Mauersberger et al. (1989}. 

(12) From trigonometric parallax of CH3OH masers, L. Moscadelli pers. comm. 

"; Source names (Col. 1); Source positions in equatorial coordinates (Cols. 2 and 3); Observed offset positions (Cols. 4 and 5); Offsets of the 
reference positions used for total power measurements (Cols. 6 and 7); Adopted distances from the Sun and the galactic centre (Cols. 8 and 9); 
Reference to these values (Col. 10). 



inner galactic disc are consistent with those of Penzias {1981), 
while the ratios from the galactic centre region are smaller. (2) 
There probably is a gradient with increasing 18 0/ 17 <9 ratios as a 
function ofRcc out to galactocentric distances of 16 kpc. (3) The 
solar system value is significantly larger than the ratios from the 
inner Galaxy, but not much higher than those from the far outer 
Galaxy. (4) An observational bias due to distance related effects 
is not apparent. 



Fig.QJ shows an increasing scatter in the 7=1-0 ratios as a 
function of Rqc- While in the inner Galaxy ratios are 5=4, there 
are a few higher ratios observed in clouds near the solar circle. 
Two (50%) of the four far outer Galaxy (FOG) clouds are also 
characterised by a high ratio. The 7=2-1 ratios show a more 
systematic increase with Rcc from about two near the galac- 
tic centre to about five at 7?gc=16 kpc. The 7=3-2 ratios may 
follow this trend, but the ratios are generally smaller. The dif- 
ferent behaviour of the 7=1-0, 2-1, and 3-2 transitions may 
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Table 3. 

Source 
(Aa,A<5) 

(arcsec) 



Observed line parameters". 



rms 
dv 



I rms 
(51 dv 



rms I rms 

dv (51 dv 

(K and Kkms- 1 ) 
(kms -1 and Kkms -1 ) 



C 18 O(l-0) C 18 0(2-l) C ls O(3-2) C 17 O(l-0) 



rms I rms I 
dv ,51 dv (51 



(1-0) (2-1) (3-2) 



(1) 


(2) 


(3) 


d) 


(5) 


(6) 


(7) 


(8) 


(9) 


WB89 380 


0.058 


0.99 


0.016 


2.89 






0.022 


0.29 


(0,0) 


0.107 


0.07 


0.107 


0.02 






0.104 


0.02 


WB89 391 


0.072 


0.58 


0.015 


1.60 






0.012 


0.078 


(0,0) 


0.107 


0.07 


0.107 


0.01 






0.104 


0.012 


DDT94 Cloud 1 


0.018 


0.20 


0.024 


0.15 






0.050 




(120,-170) 


0.107 


0.01 


0.107 


0.01 






0.104 




W 3 


0.058 


3.26 


0.075 


7.07 










(0,0) 


0.107 


0.07 


0.107 


0.09 










W 30H 


0.064 


11.8 


0.041 


20.3 


0.051 


10.2 


0.041 


3.19 


(0,0) 


0.107 


0.10 


0.107 


0.08 


0.284 


0.14 


0.104 


0.07 


WB89 437 


0.049 


1.68 


0.015 


2.86 






0.019 


0.30 


(0,0) 


0.107 


0.06 


0.107 


0.02 






0.104 


0.02 


DDT94 Cloud 2 


0.016 


0.14 


0.024 


0.20 






0.028 




(50,35) 


0.107 


0.01 


0.107 


0.01 






0.104 




WB89 501 


0.043 


0.57 


0.016 


1.83 


0.130 


2.37 


0.022 


0.17 


(0,0) 


0.107 


0.05 


0.107 


0.02 


0.142 


0.17 


0.104 


0.02 


Ori KL 


0.120 


13.8 


0.159 


47.1 


0.210 


127 


0.117 


2.52 


(0,0) 


0.107 


0.23 


0.107 


0.31 


0.285 


0.7 


0.104 


0.15 


OriKL 


0.242 


5.16 


0.098 


11.6 


0.186 


16.9 


0.090 


1.70 


(30,-240) 


0.107 


0.25 


0.107 


0.10 


0.285 


0.31 


0.104 


0.09 


Ori KL 


0.237 


10.4 


0.252 


19.3 


0.207 


35.3 


0.146 


3.25 


(30,30) 


0.107 


0.25 


0.107 


0.26 


0.285 


0.35 


0.104 


0.15 


OriKL 


0.239 


3.54 


0.074 


10.5 


0.200 


18.4 


0.071 


1.20 


(0,-150) 


0.107 


0.28 


0.107 


0.09 


0.285 


0.38 


0.104 


0.08 


NGC 2024 


0.134 


7.45 


0.061 


13.1 


0.222 


14.2 


0.080 


2.06 


(0,0) 


0.107 


0.13 


0.107 


0.06 


0.285 


0.37 


0.104 


0.08 


L 134N 






0.038 


1.29 


0.078 


0.340 


0.029 


0.46 


(-84,-48) 






0.053 


0.02 


0.142 


0.038 


0.052 


0.02 


Sgr B2(M) 


0.160 


96.7 


0.196 


136.3 


0.101 


263.2 


0.089 


33.1 


(0,0) 


0.107 


0.47 


0.107 


0.58 


0.569 


0.7 


0.208 


0.26 


Sgr B2(M) 


0.118 


64.9 


0.163 


72. (, 


0.102 


192.3 


0.071 


23.2 


(0,15) 


0.107 


0.33 


0.107 


0.44 


0.569 


0.6 


0.208 


0.19 


W 33 






0.229 


46.1 


0.087 


76.0 


0.257 


11.5 


(0,0) 






0.107 


0.34 


0.285 


0.21 


0.104 


0.38 


W 33 


0.111 


50.4 


0.093 


67.6 


0.086 


91.6 


0.063 


14.5 


(0,24) 


0.107 


0.16 


0.107 


0.14 


0.285 


0.21 


0.104 


0.09 


G34.3+0.2 


0.090 


40.2 


0.144 


63.4 


0.061 


101.2 


0.169 


10.8 


(0,0) 


0.107 


0.13 


0.107 


0.15 


0.569 


0.2 


0.104 


0.20 


W 49 


0.154 


38.9 


0.175 


90.2 


0.077 


163.1 


0.079 


11.0 


(0,0) 


0.107 


0.35 


0.107 


0.39 


0.569 


0.4 


0.208 


0.18 


W 51d-a c 


0.153 


11.5 


0.068 


12.8 


0.080 


9.64 


0.092 


3.12 


(0,0) 


0.107 


0.15 


0.107 


0.07 


0.285 


0.14 


0.104 


0.09 


W 51d-b d 


0.153 


21.0 


0.068 


34.0 


0.080 


45.2 


0.092 


6.18 


(0,0) 


0.107 


0.22 


0.107 


0.11 


0.285 


0.19 


0.104 


0.13 


DR 21 


0.099 


3.27 


0.083 


9.45 






0.049 


1.05 


(0,0) 


0.107 


0.15 


0.107 


0.13 






0.104 


0.07 


DR. 21cal 


0.112 


12.0 


0.120 


25.2 


0.098 


10.2 


0.045 


3.51 


(0,0) 


0.107 


0.17 


0.107 


0.18 


0.285 


0.25 


0.104 


0.07 


DR. 210H 


0.106 


14.1 


0.101 


24.2 


0.102 


30.5 


0.040 


3.12 


(0,0) 


0.107 


0.17 


0.107 


0.15 


0.285 


0.20 


0.104 


0.06 


NGC 7538 


0.040 


8.09 


0.050 


18.7 


0.042 


37.4 


0.030 


2.88 


(0,0) 


0.107 


0.07 


0.107 


0.08 


0.285 


0.11 


0.104 


0.05 



C 17 0(2-l) 
(10) (11) 



0.019 
0.104 
0.016 
0.104 
0.040 
0.104 
0.076 
0.104 
0.034 
0.104 
0.013 
0.104 
0.019 
0.104 
0.019 
0.104 
0.160 
0.104 
0.083 
0.104 
0.117 
0.104 
0.055 
0.104 
0.054 
0.104 
0.025 
0.052 
0.166 
0.104 
0.114 
0.104 
0.173 
0.104 
0.052 
0.104 
0.176 
0.104 
0.092 
0.104 
0.063 
0.104 
0.063 
0.104 
0.039 
0.104 
0.045 
0.104 
0.038 
0.104 
0.030 
0.104 



0.57 
0.02 
0.31 
0.01 



2.02 
0.09 
6.06 
0.05 
0.69 
0.02 



0.38 
0.02 
8.58 
0.18 
3.10 
0.08 
5.96 
0.12 
2.88 
0.07 
3.55 
0.06 
0.44 
0.02 
58.1 
0.49 
36.6 
0.30 
14.8 
0.25 
24.4 
0.07 
26.1 
0.26 
22.0 
0.21 
4.11 
0.06 
9.86 
0.09 
2.33 
0.07 
7.86 
0.07 
7.91 
0.06 
4.79 
0.05 



C 17 0(3-2) 
(12) (13) 



0.052 2.31 

0.556 0.15 

0.058 14.7 

0.279 0.15 



0.039 
0.139 
0.099 
0.278 
0.090 
0.278 
0.096 
0.278 
0.093 
0.278 
0.092 
0.278 
0.015 
0.139 
0.074 
0.556 
0.079 
0.556 
0.061 
0.278 
0.064 
0.278 
0.049 
0.556 
0.064 
0.278 
0.063 
0.278 
0.063 
0.278 



0.054 
0.278 
0.051 
0.278 
0.082 
0.278 



53.5 
0.28 
5.07 
0.15 
12.0 
0.16 
5.92 
0.16 
4.09 
0.15 
0.095 
0.006 
136.7 
0.5 
69.0 
0.49 
32.8 
0.15 
40.3 
0.02 
41.1 
0.14 
53.6 
0.21 
2.96 
0.11 
16.5 
0.15 



15.0 
0.13 
10.0 
0.14 
11.3 
0.22 



(11) 
3.57 
0.36 
7.78 
1.52 
>1.75 
(3a) 



3.87 
0.09 
5.86 
0.38 
>2.22 
(3a) 
3.51 
0.52 
5.73 
0.35 
3.18 
0.23 
3.35 
0.17 
3.09 
0.32 
3.79 
0.16 

3.71 t 
0.27 
3.06 
0.03 
2.93 
0.03 

3.98 k 
0.17 
3.64 
0.02 
3.89 
0.07 
3.70 
0.07 
3.86 
0.13 
3.56 
0.08 
3.26 
0.26 
3.57 
0.08 
4.73 
0.10 
2.94 
0.06 



18 0/C 17 ( 
(15) 
5.31 
0.19 
5.40 
0.18 

>2.08 
(3c) 
3.00 
0.17 
3.51 
0.03 
4.34 
0.14 

>6.06 
(3a) 
5.04 
0.27 
5.75 
0.13 
3.92 
0.10 
3.39 
0.09 
3.82 
0.08 
3.86 
0.00 
3.07 
0.15 
2.46 
0.02 
2.07 
0.02 
3.20 
0.06 
2.!)') 
0.01 
2.54 
0.03 
4.29 
0.04 
3.26 
0.05 
3.61 
0.03 
4.25 
0.14 
3.30 
0.04 
3.20 
0.03 
3.91 
0.04 



(16) 



2.86 
0.03 



>5.95 
(3a) 
2.49 
0.02 
3.49 
0.12 
2.89 
0.05 
3.25 
0.10 
3.63 
0.16 
3.75 
0.48 
2.02 
0.01 
2.92 
0.02 
2.43 
0.01 
2.38 
0.01 
2.58 
0.01 
3.19 
0.02 
3.41 
0.14 
2.87 
0.03 



2.81 
0.03 
3.19 
0.05 
4.09 
0.05 



rms I 
dv (SI 



13 CO(1-0) 



rms I rms I 

dv ,51 dv (51 

(K and Kkms -1 ) 
(kms -1 and Kkms -1 ) 

~ 3 CO(2-l) C 18 0(1 0) 



I (1 0) (2 1) 



C ls O(2-l) 13 CO/C 18 



(17) 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 


(26) 


0.021 


14.7 


0.151 


20.8 


0.024 


1.27 


0.022 


2.79 


11.6 


7.46 


0.106 


0.02 


0.136 


0.19 


0.107 


0.03 


0.107 


0.02 


0.3 


0.09 


0.023 


7.63 


0.124 


10.3 


0.018 


0.61 


0.018 


1.45 


12.5 


7.10 


0.106 


0.02 


0.136 


0.13 


0.107 


0.02 


0.107 


0.02 


0.2 


0.13 


0.054 


54.3 


0.174 


86.7 


0.065 


3.76 


0.118 


9.18 


14.4 


9.44 


0.272 


0.12 


0.425 


0.47 


0.273 


0.12 


0.427 


0.29 


0.5 


0.30 


0.051 


75.8 


0.081 


85.0 


0.048 


11.8 


0.071 


27.2 


6.42 


3.13 


0.272 


0.13 


0.425 


0.29 


0.273 


0.13 


0.427 


0.23 


0.07 


0.03 


0.028 


15.2 


0.079 


15.4 


0.023 


1.44 


0.026 


3.12 


10.6 


4.94 


0.106 


0.04 


0.136 


0.11 


0.107 


0.03 


0.107 


0.03 


0.2 


0.05 


0.024 


11.5 


0.182 


10.0 


0.022 


0.75 


0.027 


1.69 


15.3 


9.47 


0.106 


0.03 


0.136 


0.23 


0.107 


0.02 


0.107 


0.03 


0.4 


0.22 


0.143 


144.7 


0.403 


598.4 


0.048 


15.5 


0.190 


70.9 


9.34 


8.44 


0.272 


0.4 


0.425 


1.8 


0.273 


0.15 


0.427 


0.72 


0.09 


0.09 


0.074 


63.4 


0.108 


96.8 


0.072 


5.52 


0.113 


13.2 


11.5 


7.33 


0.272 


0.12 


0.425 


0.23 


0.273 


0.09 


0.427 


0.23 


0.2 


0.13 


0.056 


79.2 


0.111 


154.9 


0.068 


11.7 


0.107 


32.5 


6.76 


4.77 


0.272 


0.09 


0.425 


0.2 


0.273 


0.11 


0.427 


0.22 


0.06 


0.03 


0.079 


40.2 


0.109 


75.1 


0.074 


4.80 


0.104 


12.7 


8.38 


5.91 


0.272 


0.13 


0.425 


0.22 


0.273 


0.14 


0.427 


0.21 


0.25 


0.10 


0.072 


55.8 


0.092 


60.4 


0.059 


6.74 


0.099 


12.0 


8.28 


5.03 


0.272 


0.072 


0.425 


0.19 


0.273 


0.09 


0.427 


0.19 


0.12 


0.08 


0.167 


6.82 


0.122 


4.19 


0.190 


1.01 


0.114 


1.40 


4.15 


2.99 


0.027 


0.06 


0.136 


0.09 


0.027 


0.07 


0.137 


0.08 


0.18 


0.19 


0.045 


487.8 


0.255 


627.9 


0.040 


111.4 


0.098 


277.6 


4.38 


2.26 


2.72 


0.7 


0.425 


1.7 


2.73 


0.6 


0.427 


0.6 


0.02 


0.01 


0.037 


399.9 


0.213 


363.7 


0.024 


77.4 


0.092 


101.0 


5.17 


3.00 


2.72 


0.4 


0.425 


1.4 


2.73 


0.33 


0.427 


0.5 


0.02 


0.02 


0.148 


163.5 


0.441 


176.5 


0.123 


43.7 


0.217 


84.1 


3.71 


2.10 


0.213 


0.3 


0.425 


1.2 


0.213 


0.21 


0.427 


0.60 


0.02 


0.02 


0.162 


169.0 


0.292 


171.1 


0.123 


48.1 


0.278 


84.9 


3.51 


2.02 


0.213 


0.3 


0.425 


0.8 


0.213 


0.21 


0.427 


0.77 


0.02 


0.02 


0.175 


118.3 


0.257 


150.9 


0.116 


42.2 


0.284 


87.0 


2.80 


1.73 


0.213 


0.4 


0.425 


0.8 


0.213 


0.24 


0.427 


0.83 


0.02 


0.02 


0.139 


240.4 


0.224 


436.8 


0.107 


38.1 


0.263 


103.6 


6.31 


4.22 


0.213 


0.5 


0.425 


1.0 


0.213 


0.35 


0.427 


1.2 


0.07 


0.05 


0.130 


70.2 


0.183 


77.6 


0.107 


10.0 


0.203 


16.5 


6.62 


4.70 


0.213 


0.19 


0.425 


0.38 


0.213 


0.10 


0.427 


0.42 


0.10 


0.12 


0.130 


153.6 


0.183 


221.7 


0.107 


21.7 


0.203 


56.4 


7.08 


3.93 


0.213 


0.3 


0.425 


0.5 


0.213 


0.2 


0.427 


0.6 


0.07 


0.04 


0.087 


46.8 


0.197 


84.6 


0.060 


4.24 


0.121 


11.9 


11.0 


7.11 


0.213 


0.18 


0.425 


0.57 


0.213 


0.12 


0.427 


0.35 


0.3 


0.21 


0.150 


97.2 


0.234 


153.1 


0.116 


14.9 


0.280 


40.6 


6.52 


3.77 


0.213 


0.35 


0.425 


0.72 


0.213 


0.27 


0.427 


0.93 


0.12 


0.09 


0.141 


80.9 


0.253 


116.7 


0.114 


14.3 


0.281 


27.5 


5.66 


4.24 


0.213 


0.29 


0.425 


0.74 


0.213 


0.24 


0.427 


0.82 


0.10 


0.13 


0.119 


103.2 


0.146 


167.1 


0.086 


11.3 


0.191 


29.3 


9.13 


5.70 


0.213 


0.3 


0.425 


0.6 


0.213 


0.25 


0.427 


0.79 


0.20 


0.16 



Each source is cliararl r-nzcd by two lincH. First line. Col. 1: source nanie. Cols. 2 to 13: rms nuise levels and initialed iul ctisii i<>s (both in T^) for tlio l hrcc ground rul alioiiHi transitions intoiisitifs (of C lfi O 
and of C !7 0). Cols. 14 to 16: derived 7(C 18 0)//(C 17 0) integrated intensity ratios for the ,7=1-0, 2-1, and 3-2 transitions. Following Penzias (1981), the ratios were corrected for the difference in frequency 
of both isotopologues (increasing the ratios by a factor 1.047). Second line: Offset position, velocity resolution used and the uncertainties in the integrated intensities and line ratios. 

Cols. 17 to 26: the same information as the Cols. 2 to 16, but for the simultaneously observed 13 CO and C 18 (1-0) and (2-1) transitions (the spectra of these C 18 observations are not shown in Fig. A.l. 
6 For L 134N and W 33 (0,0) using May 2000 C 18 O(l-0) and August 1995 C 17 0(l-0) spectra; c : for the 45-55 kms" 1 component of W 51d; d : for the 55-75 kms" 1 component of W 51d 
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Wouterloot et al.: Interstellar 18 0/ 17 ratios 



J=l-0 



J=l-0 



o 

CO 




10 

Rgc (kpc) 

Fig. 1. C 18 0/C 17 abundance ratios, corrected for the difference 
between the C I8 and C n O frequencies but not accounting for 
optical depth effects, in the a) 7 — 1-0, b) 2-1, and c) 3-2 tran- 
sition as a function of distance from the galactic centre (Rgc)- 
Crosses show our new results. Open triangles indicate 3<x lower 
limits for DDT94 Cloud 1 and 2 that motivated this study. Open 
circles represent 7=1-0 line results from Penzias ([1981) in a) 
and 7=2-1 line results from Hofner et al. (2000) in b). Error 
bars (see Table [3) have been omitted to avoid confusion. The 
solar system value of 5.5 is marked by an asterisk. 



of diffuse gas in the high pressure environment of the galactic 
centre region, while sources further out may show large amounts 
of fractionated diffuse gas with enhanced 13 CO abundances (e.g., 
Watson et al. [19761 Bally & Langer [19821 . 

To directly determine optical depths, we fitted the relative 
strengths of the two C n O (1-0) groups of hyperfine components 
for sources with sufficiently small linewidths. In most cases, i.e. 
for WB89 391, WB89 437, WB89 501, Ori-KL (30,-240), and 
L134N, the optical depth of the C 17 O(l-0) line is <0.1. Higher 
values for the C 17 O(l-0) optical depth were found towards 
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Fig. 2. The same as Fig.Q] but with distance from the Sun plotted 
instead of ^gc- 



WB89 380 (0.78+0.72), NGC 2024 (0.82+0.38; this source 
shows more than one velocity component, but a two-component 
fit did not provide a reliable value), Ori-KL (30,30) (0.76+0.34), 
and Ori-KL(0,-150) (0.27±0.53). In all of these cases, uncer- 
tainties are high. A one velocity component fit to the C 17 0(2-l) 
transition could only be obtained towards L134N: t=0.53±0.1 1. 
Assuming r[C 18 O(l-0)] ~ 4xt[C 17 O(1-0)], these results sug- 
gest that t(C 18 O(1-0)) is commonly but not always smaller than 
unity and that t[C 1s O(2-1)] > t[C 18 O(1-0)]. Figs.|3p supports 
this finding, indicating no saturation effects in the 7=1-0 lines. 
C ls O/C 17 appears to be independent of 13 CO/C 18 0. The op- 
posite appears, however, to hold for the 2-1 lines. Here small 
13 CO/C 18 ratios are accompanied by small C 18 0/C n O values, 
suggesting some degree of 13 CO and C ls O 7=2-1 line satura- 
tion. 

Towards all sources we observed one or a few positions only. 
A map of the Orion-KL cloud was, however, made by White & 
Sandell ( [19951 ) in the C 18 and C n O (2-1) transitions. Although 
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Fig. 3. B CO/C 18 integrated intensity ratios in the a) 7 = 1-0 
and b) 2-1 transition as a function of distance from the galactic 
centre (Rgc)- Open squares indicate our new results. The values 
for the solar system of 5.5 (circle with central dot) and for the 
solar neighbourhood of 8.0 (asterisk) are also indicated, c) The 
C 18 0/C 17 integrated intensity ratio as a function of 13 CO/C 18 
for the 7=1-0 and d) J =2-1 transition. 



their primary goal was not the determination of 18 0/ 17 abun- 
dance ratios, their Fig. 8 might indicate a surprising trend to- 
wards larger C 18 0/C 17 ratios at larger extinctions or C 18 col- 
umn densities that deserves further study. Our value for Ori- 
KL(0,0) (both for the 7=1-0 and 2-1 lines) is consistent with 
the ratios found by White & Sandell d 19951 1 at high A v . 

4.2. LVG calculations 

4.2.1 . Models with C 18 0/H 2 = 1 .7x1 0" 7 

For all sources the combined data sets (excluding 13 CO, because 
some of it probably originates from a more diffuse lower den- 
sity environment) have been analysed with an LVG code, in the 
same way as done by Wouterloot et al. d2005t for the p Oph 
cloud. We used rates by Flower (2001) for collisions with H2 
and assumed an ortho/para FL abundance ratio of three (our re- 
sults are not sensitively dependent on this value). For L134N, a 
cold dark cloud, the same parameter space was used as for p Oph 
in Wouterloot et al. (f2005l > (5 K < 7 kin < 35 K, 10 3 cirT 3 < n(H 2 ) 
< 10 6 cm 4 , and 2.0 < C 18 0/C 17 < 6.0). For the other sources 
7\i n ranged from 5 K to 1 00 K or even 1 30 K (for DR2 1 cal) using 
the same range of densities and C 18 0/C 17 ratios as for L134N. 
In all cases the adopted velocity gradient was 5.0kms _I pc _1 that 
roughly reflects cloud size and measured linewidths. 

The models were compared with the observations by cal- 
culating x 1 values from the C ls O peak temperatures and the 
C I8 0/C 17 ratios, using for the latter integrated line temper- 
atures. For this procedure, the observed peak T* A temperatures 
of C I8 were converted to 7 m b temperatures except for L134N, 
where the emission is much less peaked at the observed posi- 
tion than for the other sources (see e.g., Pratap et al. 2000). The 
7=2-1 and 3-2 peak temperatures used are those of spectra that 
were convolved to the 7=1-0 resolution using results from the 
C 18 (2-1) HERA maps. For the uncertainties we assumed a cal- 




C 1S O/c"0 ratio 



Fig. 4. Results of LVG simulations modeling our sample of 
sources:^ 2 as a function of assumed C I8 0/C 17 ratio. The mini- 
mum^ 2 values are indicated by filled triangles (SgrB2), squares 
(3kpc<7?cc<10kpc), and circles (7?cc>16kpc). The calculations 
were made with steps of 0.125 (near the minimum^ 2 ) to 0.50 in 
the C 18 0/C 17 ratio. The C 18 0/H2 abundance ratio was assumed 
to be 1.7 10~ 7 (Frerking et al. 119821) . The results are similar to 
those also accounting for a galactic C 18 0/H2 abundance gradi- 
ent. 



ibration error of +10% in the line temperatures and of +14% in 
the line ratios (see Sect. 2.3). In Fig. f4j x 2 1S shown as a func- 
tion of the adopted C 18 0/C 17 ratio. This ratio was varied in 
steps of 0.125 (near the minimum^ 2 ) to 0.5. For these models 
we assumed the Frerking et al. (|T982) C 18 0/Fi2 abundance ratio 
of 1.7 10~ 7 . All sources show a well-defined minimum. The pa- 
rameters of the model with the lowest x 2 are listed in Table [4] 
Col. 4. 

The average C 18 0/C 17 ratio is 4.18+0.15. Excluding 
the two positions towards Sgr B2, which have low ratios 
(2.88+0.1 1), and the objects at /?cc>16kpc, which have high val- 
ues (5.02±0.45), the average ratio becomes 4.15+0.10. This is 
close to the values (4.11+0.14 and 4.13+0.13) derived for p Oph 
and NGC 6334 FIR II by Wouterloot et al. d2"0"05l l and Zhang et 
al. d2007), respectively. 

4.2.2. Introducing a C 18 abundance gradient 

From optical and FIR observations it has been shown that many 
metals show a considerable abundance gradient. Abundances 
decrease from the galactic centre to larger Rgc (e.g., Rudolph 
et al. [2006b . Gradients amount to -0.078 dex/kpc for 14 N/H, 
-0.051 dex/kpc for ie O/H, and -0.044 dex/kpc for 32 S/H (these 
are average values from optical and FIR measurements which 
yield slightly different results). For 12 C/H, there are only few 
data. Rolleston et al. ( 120001 ) obtained a slope of -0.070 dex/kpc 
from 80 stars in 19 open clusters (6</?GC^18kpc). Esteban et al. 
(2005) observed recombination lines towards 8 Hn regions (6.3 
< Rgc < 10.4kpc) and derived a slope of -0. 103 dex/kpc. Model 
calculations by Matteucci & Francois ( 1989) suggest a gradient 
of -0.066 dex/kpc for 12 C/H. In all cases, where measurements 
extend to large galactocentric radii, there is no sign for a change 
in the radial gradient. 

In addition to the gradients in the 12 C and 16 abundances, 
there is also a gradient in the 16 0/ 18 ratio: Wilson & Rood 
d 19941 ) give a ratio of ~250 for the galactic centre region and 



8 Wouterloot et al.: Interstellar 18 0/ 17 ratios 



Table 4. Results of LVG model calculations". 



Nr 


Source 


Offset 


c 18 o/c' 7 o 


7V„ 

1 km 




Y 1 
A 


C 18 0/C 17 


logrc(H 2 ) 


y 2 

A 


-*^GC 












(cm' 3 ) 
Win ) 






(K) 


crrr 3 




kpc 








c 


8 0/Hi 






c 


s O/H 2 Frerking+gradient 






1 


WB89 380 


0,0 


4.13 


77.5 


4.0 


3.6 


4.13 


45.0 


4.8 


3.4 


16.8 


2 


WB89 391 


0.0 


6.31 


70.0 


4.1 


3.1 


6.38 


45.0 


4.9 


2.6 


16.9 


3 


W 3 


0.0 


3.75 


25.0 


4.3 


3.0 


3.88 


25.0 


4.5 


6.2 


10.0 


4 


W 30H 


0,0 


3.88 


30.0 


4.8 


0.9 


3.88 


27.5 


4.9 


1.6 


10.0 


5 


WB89 437 


0.0 


5.13 


30.0 


4.0 


2.0 


5.13 


20.0 


4.7 


2.4 


16.2 


6 


WB89 501 


0.0 


4.5 


60.0 


3.9 


7.2 


4.5 


20.0 


4.7 


16.4 


16.3 


7 


Ori KL 


0,0 


3.88 


52.5 


4.9 


19.5 


3.88 


62.5 


5.0 


19.1 


8.84 


8 


Ori KL 


30,30 


4.0 


45.0 


5.3 


1.2 


4.13 


42.5 


5.3 


2.1 


8.84 


9 


Ori KL 


0,-150 


3.75 


45.0 


4.9 


2.4 


3.75 


52.5 


5.0 


2.4 


8.84 


10 


Ori KL 


30,-240 


4.0 


40.0 


5.0 


3.5 


4.0 


37.5 


5.0 


4.8 


8.84 


1 1 


NGC 2024 


0,0 


4.63 


20.0 


4.8 


3.8 


4.5 


20.0 


4.8 


5.7 


8.86 


12 


L 134N 


0.0 


5.0 


7.0 


4.7 


5.1 


5.0 


7.0 


4.6 


5.6 


8.41 


13 


Sgr B2M 


0,0 


3.0 


37.5 


5.1 


6.1 


3.0 


42.5 


4.3 


5.3 


0.1 


14 


Sgr B2M 


0,15 


2.75 


35.0 


4.8 


18.6 


2.75 


57.5 


3.9 


21.5 


0.1 


15 


W33 


0,0 


4.0 


30.0 


5.1 


1.8 


4.25 


30.0 


4.8 


1.1 


4.8 


16 


W33 


0,24 


5.13 


30.0 


5.5 


1.5 


4. / J 


32.5 


5.0 


1.6 


4.8 


17 


G34.3+0.2 


0,0 


4.13 


35.0 


5.3 


2.5 


4.13 


35.0 


5.0 


2.2 


5.7 


18 


W49N 


0.0 


4.13 


50.0 


5.0 


3.0 


4.25 


45.0 


4.9 


3.2 


7.8 


19 


W51d-a fc 


0.0 


4.44 


12.5 


4.7 


1.5 


4.5 


12.5 


4.4 


1.3 


6.6 


20 


W51d-r/ 


0.0 


3.88 


25.0 


4.8 


3.4 


4.0 


25.0 


4.7 


3.6 


6.6 


21 


DR21 


0.0 


3.75 


107.5 


4.7 


1.7 


3.75 


107.5 


4.7 


1.7 


8.5 


22 


DR 2 leal 


0.0 


4.0 


30.0 


5.0 


3.7 


4.0 


30.0 


5.0 


3.7 


8.5 


23 


DR 210H 


0.0 


4.63 


22.5 


4.9 


2.5 


4.63 


22.5 


4.9 


2.5 


8.5 


24 


NGC 7538 


0.0 


3.69 


40.0 


4.8 


6.8 


3.63 


40.0 


4.9 


6.7 


9.8 



Cols. 2 and 3: source and offset; Cols. 4 to 7: derived C 18 0/C 17 ratio, kinetic temperature, H 2 density and minimum x 1 ■ 
Cols. 8 to 11: the same parameters as Cols. 4 to 7, but for the case of a galactic C ls O/H2 abundance gradient (see Sects. 4.1 and 4.2). 
: W51d-a: 50kms~' component; W51d-b: 60kms~' component (see Fig. IA.lt . 



a fitted ratio of 16 0/ 18 ~ (59±12)xfi G c + (37±83) for larger 
Rcc values (2.8 < R C c < 9 kpc). 

Adopting the above mentioned 16 0/H abundance gradient of 
-0.05 1 dex/kpc, the oxygen abundance becomes a factor of 2.7 
higher near the galactic centre and a factor of 2.4 lower at Rcc 
= 16 kpc than in the solar neighbourhood. Similarly, the 16 0/ 18 
ratio is then factor of 2.1 lower near the centre (using a mini- 
mum 16 0/ 18 ratio of 250) and a factor of 1.8 higher at R cc = 
16 kpc than the fitted value near the Sun, 537 (Wilson & Rood 
1994) which can be compared with the solar system value of 490 
(Anders & Grevesse \T9E% . 

Applying both gradients to estimate ls O/H, the ratios be- 
come a factor of 5.7 higher at the centre and a factor of 4.4 
lower at Rcc - 16 kpc than in the solar neighbourhood. Taking 
for 12 C the mean of the gradients determined by Rolleston et 
al. (120001 and Esteban et al. d20"05l -0.0865 dex/kpc, we obtain 
12 C/H abundances that vary between 5.4 higher and 4.5 times 
lower within the same range of galactocentric radii. Since the 
gradients are uncertain and the true abundances also depend on 
chemical reactions in the clouds, the 18 0/H and 12 C/H gradi- 
ents are identical within the error limits. In our LVG calculations 
we have therefore used the ls O/H gradient for all sources to ob- 
tain a second estimate of 18 0/ 17 ratios. For large Rgc, the ex- 
trapolated lower abundances are consistent with chemical model 
calculations (Ruffle et al. 120071 1 simulating multi-line measure- 
ments of DDT94 Cloud 2. 

The results are given in Table [4] (Cols. 7 to 10) and Fig. [5] 
There is little change in the C 18 0/C 17 ratios compared to those 
without considering the abundance gradient (Col. 4 of Table 0). 
The average C 18 0/C 17 ratio becomes 4.20±0.14 instead of 
4.18±0.15. Excluding the two positions in Sgr B2, which have 
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Fig. 5. Results of LVG calculations: The C 18 0/C 17 ratio as 
a function of Rgc- a) The C ls O/H2 ratio was assumed to be 
1.7 10~ 7 (Frerking et al. 1982), or b) corrected for the galactic 
abundance gradients (see Sect. 4.2.2). 



low ratios (2.88+0.11), and the objects at /?GC>16kpc, which 
have higher values (5.03+0.46), the average ratio becomes 
4. 16+0.09. This is again close to the value derived by Wouterloot 
et al. (120051 ) forp Oph (see also Zhang et al. |2007l for NGC6334 
FIR II). There are, however, differences in the derived kinetic 
temperatures and densities. In the first set of models (with fixed 
C ls O/H2 abundance ratios) the derived densities are significantly 
smaller at 7?cc>16kpc than for 3</?GC<10kpc. This difference 
is greatly reduced in the second set of models accounting for 
C 18 0/H2 gradients (Pig. |6j>. If we assume that the H2 densities 
in star forming clumps such as those in our selected sample are 
similar at all locations in the Galaxy, this indicates that the mod- 
els with the abundance gradient are more realistic. 
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Fig. 6. Results of LVG calculations: The crosses indicate the best 
fit 7\i n and n(H2) using the Frerking et al. ( 1982 ) abundance. The 
arrows point to the results for the model using an abundance 
gradient. The arrows for the clouds at R>16kpc are dotted, those 
for Sgr B2 are dashed. Numbers indicate the sequential line-of- 
sight number, indicated in Table [4] 



4.2.3. Implications 

The two sets of models, either ignoring or accounting for ra- 
dial C ls O/H2 abundance gradients across the Galaxy, suggest the 
existence of a C 18 0/C 17 gradient from values slightly smaller 
than three near the galactic centre to values around five at large 
galactocentric radii (Fig. |5). The 18 0/ 17 ratio of Sgr B2 de- 
duced from HCO + , 3.1+0.6 (Guelin et al. 119821 ), is consistent 
with this result. Our results do suggest that the C 18 0/C 17 ratio 
at 4 to 1 1 kpc from the galactic centre is larger than the value de- 
rived by Penzias dl98U 3.65+0.15). The reason is that the 7=1- 
results from Penzias were not corrected for radiative transfer 
effects. 

18 is synthesised by helium burning on 14 N in massive stars. 
14 N is a mostly secondary nucleus so that its abundance is highly 
metallicity-dependent, with low N/C and N/O abundance ratios 
at low metallicities (e.g., Wheeler et al. 119891 . 18 may follow 
this trend. n O is mainly a product of CNO-hydrogen burning 
and is also released by stars of intermediate mass (for yields, 
see, e.g., Prantzos et al.[T996). 

To qualitatively reproduce the 18 0/ 17 ratios observed in in- 
terstellar clouds and stellar objects, we should note that in in- 
dividual dust grains from massive stars ls O overabundances of 
up to two orders of magnitude can be reached (e.g., Amari et al. 
1995 ). Low ratios, sometimes below unity, are obtained towards 
late-type stars of intermediate mass (c.f., Sect. 1). The low val- 
ues indicate that 18 is destroyed in such stars (see also Henkel 
& Mauersberger UHH Henkel et al. [T994b . 

For the Large Magellanic Cloud (LMC ), the C 18 0/C 17 Q (2- 
1) ratio has been studied by Heikkila et al. ( 1998). These Clouds 
have even lower metallicities than the far outer Galaxy and are 
characterised by C 18 0/C n O ~ 1.6, which is significantly lower 
than all measured interstellar values in the Galaxy. The LMC, 
which undergoes a phase of enhanced massive star formation, 
is metal-poor so that its 18 abundance should be low. Very 
high C 18 0/C 17 ratios are instead found in nuclear starbursts 
(e.g., Harrison et al. |1999t Wang et al. 120041 ). This may be a 
consequence of (1) high metallicities and (2) large amounts of 



ejecta from massive stars that must have enriched the interstellar 
medium. 

In the central region of our Galaxy, metallicities are also 
high. However, here CNO burning dominates resulting isotope 
ratios of C, N, and O. Apparently, a nuclear starburst with large 
numbers of massive stars has not contaminated the galactic cen- 
tre region since a long time. Products from helium burning are 
therefore underepresented in the C, N, and O isotope abun- 
dances, yielding a very low 12 C/ 13 C ratio (~25), an extremely 
large 14 N/ 15 N ratio (of order 1000) (e.g., Wilson & Rood [19941 
and the low 18 0/ 17 ratio (SS 3) proposed here. 

What remains, at first sight, puzzling is (1) that there exists a 
galactic disc 18 0/ 17 gradient and (2) that the high 18 0/ 17 ra- 
tios in the metal-poor outer Galaxy do not match the low values 
in the even more metal-poor LMC. Does this latter discrepancy 
imply that the metallicities of the outer Galaxy still permit effi- 
cient high-mass star production of 18 0? While the LMC is too 
metal-poor to make the process efficient? Or that, in spite of the 
presently visible large numbers of massive stars the LMC star- 
burst is still too young to enrich the interstellar medium with 
18 0-rich ejecta? Both possibilities appear to be farfetched. The 
first, because the metallicities are not that different (e.g., Hunter 
et al. 120071 Ruffle et al.|2007 ) to justify an extreme variation (by 
a factor of three) between measured 18 0/ 17 ratios, the second, 
because the starburst in the LMC started as early as 50xl0 6 yr 
ago (e.g., Westerlund 1990]). 

The galactic gradient and the difference between outer 
Galaxy and LMC may instead be interpreted in terms of cur- 
rent models of galacto-chemical evolution. The most accepted 
mechanism to explain the existence of abundance gradients is 
the so-called "biased infall" (e.g., Chiappini & Mateucci 1999) 
with the galactic disc being slowly formed from inside out. In 
such a scenario, 18 abundances may be low in both the outer 
disc and LMC and the difference is mainly caused by 17 0. In 
the LMC, where star formation is ongoing since at least 10 ll) yr 
(e.g., Hodge [T989), there was sufficient time to form n O pre- 
dominantly from stars of intermediate mass. The outer Galaxy, 
however, may be too young to build up a similar n O abundance, 
with timescales (and abundances) rapidly decreasing the farther 
out the measured molecular cloud is located. 

To summarise, we may have found a fundamental differ- 
ence between the metal-poor outer regions of the Galaxy and 
the metal-poor LMC. Generalizing this result, the difference im- 
plies that the abundances of the outer regions of spirals cannot be 
considered to be intermediate between those found near the so- 
lar circle and those determined in small metal-poor dwarf galax- 
ies. As a cautionary note, however, we have to emphasise that 
the 18 0/ 17 (9 ratios of the innermost and outermost galactic star- 
forming regions are still based on a very small number of targets. 
Additional data to confirm or to reject the trend of rising 18 0/ 17 
ratios with galactocentric radius would thus be highly desirable. 



5. Conclusions 

To measure interstellar 18 0/ 17 ratios, we observed three C 18 
and C n O transitions towards 25 positions in 18 galactic sources 
located at galactocentric distances Okpc < Rqc < 16.9 kpc. 
These measurements are complemented by 13 CO observations 
in the two ground rotational transitions. C 18 0/C n O line inten- 
sity ratios from the 7=2-1 and 3-2 transitions are lower than 
those from the 1-0 lines, indicating that optical depth effects 
have to be considered. This is also suggested by a comparison 
with the 13 CO (1-0) and (2-1) spectra. 



Wouterloot et al.: Interstellar 18 0/ 17 ratios 



10 

Therefore, all C ls O and C 17 observations were combined 
by means of LVG calculations. This was done assuming (1) that 
the C ls O/H2 abundance ratio is constant throughout the Galaxy 
and (2) that there is a radial C 18 0/H2 abundance gradient. These 
assumptions do not affect significantly the final C 18 0/C 17 val- 
ues, which should provide a good approximation to the 18 0/ 17 
isotope ratio. For the central region, 18 0/ 17 = 2.88+0.11. For 
^gc = 4-llkpc, the ratio becomes 4.16+0.09. In the outer 
Galaxy (16kpc < R cc < 17kpc), we find 5.03+0.46. 

The low ratio in the galactic centre region is consistent with 
a CNO-hydrogen burning dominated nucleosynthesis that is also 
characterizing the carbon and nitrogen isotope ratios. It supports 
the view that ls O is predominantly synthesised in high-mass 
stars, while 17 is predominantly a product of stars of lower 
mass. The ratio between 4 and 1 1 kpc is consistent with recent 
results from a few "local" individual clouds. It is smaller than the 
solar system ratio (5.5), suggesting that the Sun was enriched by 
material from massive stars during its formation. The 18 0/ 17 
values measured for the outer Galaxy are more difficult to inter- 
pret. They are not low, as in the case of the metal-poor massive 
star forming Large Magellanic Cloud but appear to be higher 
than anywhere else in the interstellar medium of the Galaxy. 
This and the galactic disc gradient may be explained by the small 
age of the outer galactic disc and may imply that the metal-poor 
outer reaches of spiral galaxies provide quite different environ- 
ments than similarly metal-poor dwarf galaxies. Nevertheless, 
both the Galactic centre and outer Galaxy 18 0/ 17 isotope ratios 
suffer from small number statistics so that more observations are 
needed to confirm or to reject the trend found in this study. 
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Appendix A: Spectra 
List of Objects 

'DDT94 Cloud 1' on page 
'DDT94 Cloud 2' on page |2 
'WB89 380' on pageH 
"WB89 391 'on pageg 
'DDT94 Cloud 1' on pagei 
'W3' on pagei 
'W 30H' on pagei 
'WB89 437' on pagei 
'DDT94 Cloud 2' on pagei 
'WB89 501' on pagei 
'Ori KL' on pagei 
'NGC 2024' on pagei 
'L 134N' on pagei 
'SgrB2(M)' on pagei 
'W33'on pagei 
'G34.3+0.2' on pagei 
'W 49N' on pagei 
'W 51' on pagei 
'DR 21' on pagei 
'DR210H' onpage[ 
'NGC 7538' on pagefc 
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Fig. A.l. Continued. 



